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INTRODUCTION 

Statistical properties of quantum many-body systems 
are determined by calculating thermal averages with a 
proper distribution function, which represents the mean 
number of particles occupying each of the single-particle 
states. For the three dimensional systems of noninteract- 
ing particles this distribution function is given by either 
the Fermi-Dirac (FD) or the Bose-Einstein (BE) func- 
tion, depending on whether the spin of the particles is 
half-integer or integer, respectively Q. 

In the systems with either strong interactions or of 
low space dimensionality (d < 3) the situation changes 
drastically. It turns out that the population of single par- 
ticle quantum level can be given by a function which is 
of neither FD, nor BE form. The reason for that may 
be twofold. First of all, in the space of two dimensions 
the symmetry of the many body wave function is not 
necessarily either even or odd |2(. When two identical 
particles are exchanged the wave function acquires an 
arbitrary phase factor e la , where a is so called statisti- 
cal angle. In other words, the proper symmetry group 
in this case is the braid group for which the irreducible 
representations are x{ a ) — e M - The particles obeying 
the fractional statistics are called anyons. On the other 
hand, one can imagine that the interactions between par- 
ticles are sufficiently strong, e.g. singular, that number 
of available effective single-particle occupancies of each 
of the levels depends on the number of particles already 
present in that state. This property can also lead to the 
departure from FD and BE statistics in arbitrary spatial 
dimensions |3j. 

There is no universal theory of statistical properties 
of the systems described above. However, there are 
some exactly solvable models of interacting particles for 
which the momentum distribution functions are modi- 
fied Q . One of the best known cases is the one dimen- 
sional Calogero-Suthcrland model for which the statisti- 



cal properties of quasiparticles change from BE to FD 
limits as the interaction increases 

Having in mind those exactly solvable examples and 
the fact that, in general, the statistics of elementary ex- 
citations (quasiparticles) in correlated systems can be 
arbitrary, we believe it is valuable to study the ther- 
mal properties of particles with nonstandard distribution 
functions. In particular, in this paper we address the 
question how thermodynamic properties evolve when the 
statistics changes from BE to FD limits. 

Our paper is organized as follows. In the first Sec- 
tion we introduce three interpolation schemes between 
BE and FD statistics. Then, in Section 2 and 3 we dis- 
cuss the high and the low temperature properties of such 
systems. In Section 5 we present general phase diagrams 
for those three systems of particles. 



ROUTES OF INTERPOLATING BETWEEN 
FERMI-DIRAC AND BOSE-EINSTEIN 
DISTRIBUTION FUNCTIONS 

For noninteracting bosons the momentum distribution 
function is 



n(ei) 
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whereas for ideal fermions we have 
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In these functions, and hereinafter, denotes the energy 
of a single particle state i, (3 = X/ksT, where T is the 
temperature, and /i is the chemical potential. Probably, 
the simplest modification of FD and BE functions can be 
carried out phenomenologically in the following manner 
(route 1) 
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where 7 is an arbitrary real number in the range (—1,1). 
The limit —1 corresponds to BE statistics and 7 = 1 
to FD one. For other allowed values of 7 we have an 
intermediate statistics. In particular, the limit 7 = 
corresponds to the classical, i.e. Boltzman, distribution 
function. Particles with arbitrary 7 statistics are called 
7 -oris. We show later that this ad hoc proposal rep- 
resents a high-temperature approximation to a general 
distribution function. 

As far as we know, it is impossible to derive the dis- 
tribution function © exactly [(J til- However, starting 
from the modified Hilbert-Fock space with creation and 
annihilation operators obeying so called g-mutator alge- 
bra, i.e. aiO^ — e la a^ai — for i ^ j, there is a way to 
recover the function J3J) approximately where 7 = cos a 
7] . It is also interesting to observe that the distribution 
function can be obtained if we suppose that the num- 
ber of ways of putting Ni particles into Gi states is given 
by 
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Then, using the standard combinatoric approach to 
quantum statistics |l| we recover eq.©- 

The thermodynamic potential S 7 (z, V, T) for 7-ons is 
given by 



& r (z,V,T) = '[[(l+<yzeP e <) 1 



h 



(5) 



where z = e^** is the fugacity. In eas. (|4l5(l we supposed 
that 7^0. The case of 7 = must be treated separately. 
The equation of state for such gas of ideal 7-ons is 



pV 



= lnS 7 (z,y,T) = i^ln(l+7ze' 3e -), (6) 



where p is the pressure and V is the volume of the system. 
In order to eliminate z from this formula we need another 
equation which determines the total number of particles 
N 



N = z-^lnS^z^T) = 



(7) 

In deriving the last equation we could also convince our- 
selves that the thermodynamic potential JSJ is really con- 
sistent with the distribution function (J3J. Eqs. JBJ) and 
describe the whole thermodynamics of 7-ons. We are 
going to explore them in the next sections. 

On the other hand, one can interpolate between FD 
and BE limits using the state counting arguments. This 
provides us with the second possibility (route 2). Namely, 
let the number of single particle state Di in the i-th 
quantum level depends on the number of particles in this 
group according to the formula [3J, |8| 



Also, suppose that the total number of states of N-body 
system for a given configuration {N} is 



w({n 1 })=h 



(d - g(Nj - 1) - 1 + Nj)\ 
(iV,)!(G l - ff (iV i -l)-l)! : 



(9) 



then using the Boltzman equation for the entropy S = 
k B hxT(e), where T(e) = J2{N-} W ({ N i}) is the total 
number of states, we can derive the momentum distri- 
bution function n fl (e,) Q. In eqs. JBJ and © Gi is 
the bare number of single particle states (when no par- 
ticles are present), g is so called statistical interaction, 
or statistical parameter, which in fact parameterize the 
distribution function. We note that for g — 0, the num- 
ber of states is the same as for the ideal bosons. The 
case 3 = 1 corresponds to free fermions. However, for 
< g < 1, we can again obtain an arbitrary distribution 
function. Physical interpretation of g is follows. Ideal 
bosons do not obey exclusion principle, so the number 
of states in the i-th energy group would not change if 
we put another particle into this system and, therefore, 
Di({Nj}) — Gi. Fermions obey Pauli exclusion princi- 
ple, which states that at most only one fermion can be 
present in a given quantum state. Then, if we put another 
fermion into the i-th group the number of accessible state 
must decrease by one. Hence, Di({Nj}) = Gi — Ni + 1 
and g = I. For arbitrary g between zero and one we 
have a partial exclusion principle and equation di- 
rectly leads to fractional exclusion statistics in arbitrary 
spatial dimensions. Explicitly, the distribution function 
for any g is given by the following transcendental equa- 
tion [3 

ntWePto-ri = (1 + {l-g)n g {e i )) 1 -o{l-gn g {e i )Y. 

( 10 ) 

Of course, taking g = or g = 1 in this equation we re- 
cover distribution functions and (J2J. Particles obey- 
ing statistics with arbitrary g are called g-ons. From the 
last equation we find that the distribution function can 
be written in the following form 



n ff (e*) 



(11) 



w g (ei)+g 
where w g (ei) must satisfy the equation 

e ^-^) = ( % (e i ) + l) 1 - ff K(ei)) S - (12) 

As in the standard statistical mechanics we can find 
the thermodynamic potential for g-ons E, g (z,V,T) and, 
hence, the equation of state is 



pV 



^ ln f'l + (l-.9K( e 



1 -gng{ti) 



D i ({N j }) = G i -g(N i -l). 



(8) 



The equation for the total number of g-ons is 

N -J2^ — ■ 

V w g (e z )+g 



(13) 



(14) 
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The function n g (ei) and w g (ci) are given by eqs. (|10fl and 
(|12|l respectively. Again, these functions i|13|) and i|14|) 
with (|1(J[I and H12|) determine the whole thermodynamics 
of g-ons. 

Finally, there is another way of changing FD and BE 
statistics (route 3). Namely, suppose that the number 
of particles in each single particle quantum state can be 
equal to p, where p is a natural number, i.e. 



Ni = 0, 1, 2, ...p. 



(15) 



The case p = 1 corresponds to free fermions and p — oo 
gives BE statistics. The thermodynamic potential for a 
given p is 



p 1 p(p+i)/9(£ ! -m) 

s^T)=n£^ (e< -' K =nV 

i rii— i 



(16) 

To derive this we used the method of grand canoni- 
cal ensemble. With the aid of the property n p (ei) — 
zd / dzln'E.p we can obtain the distribution function 



n p {ei) 



1 



p+l 



eP{€i-n) - i e (p+i)^(e 4 -/i) _ \ ' 
The equation of state for a gas of so called p-ons is 



pV 



hi- 



1 



r P+l p -f3U(p+l) 



k B T l-ze-Pe* 
and equation for their number is 



N 



E 



i 



p+i 



oP(ti-n) _ i e (p+i) j 9(e < -n) _ i 



(17) 



(18) 



(19) 



Let us summarize this Section. We introduced three 
schemes of modified quantum statistics and, as a result, 
we obtained three families of new nonstandard statistical 
distribution functions which are parameterized by some 
numbers (-y,g,p). In these cases the statistics interpolate 
between FD and BE limits as the statistical parameters 
are changed. In the next sections we are going to discuss 
thermal properties of such systems as a function of these 
parameters. 

HIGH TEMPERATURE PROPERTIES 

Let us consider first the ideal 7-ons with parabolic dis- 
persion relation tk — moving in the d-dimensional 
space. Then, we can expand the functions © and JZJ in 
powers of 72, and performing the simple Gaussian inte- 
grals we find that 



N 
V 



1 1 



(20) 



(21) 



where the function fk(l,z) is defined as 

/*(7,*)=E(-l) i+1 T ( 22 ) 



A = 2t: h 2 /rnksT is the thermal wavelength, and n is 
the particle density. Of course, one can simply check 
that these equations reduce to those for ideal fermions or 
bosons as 7 = 1 or 7 = — 1 respectively. 

Now, considering the high temperature limit, i.e. a 
limit for which nX d j C 1, we can solve eq. © finding 
the fugacity 



z ~ n\ a + 



d , 7 («A rf ) 2 



2 d/2 



(23) 



Then, substituting this into ea. (|20|l we obtain the high- 
temperature (low-density) limit of the equation of state 



knT 



= n(l + A1n-\ ), 



where is a second virial coefficient 



X" 



2 d/2+l ' 



(24) 



(25) 



The character of the first quantum correction to the clas- 
sical gas depends on the sign of 7. For 7 > 0, is 
positive and, effectively, particles repel each other. In 
this case the pressure is higher than for the classical gas. 
For 7 < the second virial coefficient is negative which 
means that the pressure is lowered. In this range of sta- 
tistical parameter 7 particles tend to collapse because 
they effectively attract each other. Roughly speaking, 
for 7 > the particles resemble fermions and for 7 < 
they look like bosons. At the boundary (7 = 0) between 
these two limits there is no quantum correction to the 
classical behavior. 

We can also calculate the internal energy U = (d/2)pV 
and, hence, the specific heat of 7-ons at high tempera- 
tures. The result is 



C v 



ksd 



n(l+^n(l-rf/2) + •••)■ 



(26) 



In the one dimensional case the specific heat is decreas- 
ing function of temperature for positive 7. For d = 3 
the behavior is opposite whereas, for d = 2 the specific 
heat does not depend on temperature. Also, we see that 
as T — > 00 the specific heat goes to its classical value 
(d/2)R, where R is a gas constant {R = nks), and this 
result is independent of 7. 

Next, we examine the high temperature properties of 
g-ons. The momentum distribution function for them 
is not given explicitly but via the solution of eq. I jlOjl . 
However, we can find exactly this solution in term of 
Taylor power expansion |l0| . To do this we solve eq. 



4 



(|10|l by utilizing the Legendre inversion theorem. The 
distribution function is in the form of the expansion 



n g (ei) 



(mg + 9 - m) m (-1) 
2— 1 ml 

m—l 



t. 



m+1 ' 



(27) 



where £j = e^ ei_ W and (a) m — a(a + 1) ■■■(a + m— 1) = 
T(a + m)/T(a). For details see Appendix A. Then it is 
straightforward to find the equation of state for ideal g- 
ons 



V 



1 



2fd/2+l(g,z), 



k B T X° 

and the equation for the total number of particles 
N 1 f i \ 

n= v= jdfd/2(g,z)- 



(28) 



(29) 



Again, A is the thermal wavelength and function fk(g, z) 
is defined as follows 



fk(g,z) 



OC 

£ 

m=0 



(-1)' 



(gm + g - m) m z 



m+l 



(m+ l) fc ' 



(30) 



Eliminating z from eq. 1)28(1 with the aid of eq. ((29(1 we 
find the equation of state in the high temperature limit 
in the form 



V 



k B T 



= n(l + A%n-\ ), 



where the second virial coefficient for g-ons is 

A'' 



(2.9 ~ 1) 



2 d/2+r 



(31) 



(32) 



For g > 1/2 this coefficient is positive, so g-ons are like 
fermions. However, for g < 1/2 they resemble bosons in 
that sense that they tend to condense. We also see that 
identifying statistical parameters 7 and g, i.e. 



7 = 2.9 - 1, 



(33) 



these two sorts of particles obeying fractional statistics 
are identical at high temperatures. Also, the whole re- 
sults for the internal energy and the specific heat for 7- 
ons could be repeated for <?-ons if we use the relation 
(ESI. 

Finally, we examine the high temperature properties 
of p-ons. Again, the distribution function l(17|) can be ex- 
panded in power series and Gaussian integrals can be per- 
formed easily for the ideal-gas case. We find the equation 
of state and equation for particle number in the forms 



1 



k B T X d 



TTi9d/2+ig{P,z), 



N 
V 



9d/2(p, 



(34) 



(35) 



In the present case the function g^ (p, z) is following 



m—l 



m(p+l) 



If we take p > 1 then the fugacity z is 



(n\ d f 

2 d/2 



(37) 



where we dropped all higher order terms. For p = 1 the 
sign in the second term on the right hand side of eq. 1(3 7|) 
would be positive. Substituting 1(37(1 into 1(34(1 and taking 
only terms of the order n 2 we have 



= n(l + A p 2 n + ■ ■ •), 



where 



A 'P - A 

A 2 - 2 d/2+l ' 



(38) 



(39) 



We see that the second virial coefficient for p-ons does not 
depend on p and is always negative for p > 1. Therefore, 
we conclude that at high temperatures these particles 
resemble ordinary bosons. 

To summarize this Section, we notice that for the 
first two statistical distribution functions we obtained the 
same behavior in the high temperature limit, whereas in 
the third case the result was different. In other words, 
<7-ons and 7-ons behave similarly whereas p-ons always 
have the boson type corrections to classical limit. In the 
next section we check the low temperatures properties of 
those systems. 



LOW TEMPERATURE PROPERTIES 

Low temperature properties of 7-ons must be discussed 
separately for positive and negative 7. Let us consider 
first the 7 > case. 

At T = the momentum distribution function reduces 
to the step function 



%(e») = -6(n - €i). 
7 



(40) 



All states below the Fermi level are occupied with the 
mean population 1/7. As 7 — » this number diverges. 
For finite 7 we see that the averaged number of particles 
at each quantum level could be fractional and greater 
than one. The states above the Fermi level are empty. 

For free particles we can find the Fermi energy — 
fi(T = 0). Namely, using eq. ((40(1 we calculate that for a 
given density of particles 



2Q 



(41) 
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where n = (2^r { d/i 



(■p-) d / 2 - The length of the Fermi The equation for total number of particles is 



vector |kp 



2ra£j 



depends on the value of 7 and so 
does the volume enclosed by the Fermi surface. In par- 
ticular we see that as 7 — > the Fermi vector disappears 
and all particles are on the lowest quantum level. 

The existence of the Fermi surface implies that the low 
temperature properties of a system are affected only by 
the low energy particle-hole excitations close to it. In 
other words, only the particles from the vicinity of the 
Fermi surface give contributions to the thermodynamics 
in this limit. Therefore, we perform the low temperature 
(Sommerfeld type) expansion for any thermal quantity. 
However, in order to do this we have to make some ex- 
tra modification in the distribution function J2J). Let us 
rewrite it as follows 



1 



1 



1 



3 /3(£-m) 7 -f -f e /3(e-p*) + 1 ' 



(42) 



where /i* = /1 + &bT In 7. Now, we can expand the ther- 
modynamic functions around the point /i*. The general 
scheme of calculation is presented in Appendix B. Here 
we quote only the final results. 

The lowest corrections to the chemical potential are 



k B T\n^ — 



2 1 



6 24 



bl 



(k B Tf 



(43) 



For ordinary fermions (7 = 1) we had the first correction 
to the chemical potential proportional to T 2 . Now, the 
linear in temperature correction has appeared. This is a 
direct consequence of the modified form of the distribu- 
tion function ©. However, this term disappears in the 
formula for the internal energy. Namely, as was shown in 
Appendix B the internal energy is 



U 



id\"l\ 

2n 



{k B Tf 



(44) 



The general form is similar to that for fermions. The only 
difference is the factor 1 /j which enhances the density of 
states at the Fermi level. Also, the specific heat 



C v = K-yT 



(45) 



(1/7) 



2/d 



is linear function of T. The coefficient K, 
diverges as 7 — > 0. 

So, we proved that the ground state and the low tem- 
peratures properties of 7-ons for 7 > resemble proper- 
ties of ordinary fermions. However, it is not the case for 
7 < 0. 

For negative statistical parameter 7 the distribution 
function J2J is singular when + /eg Tin 7 — /i = 0. This 
may be a signal of the Bose-Einstien condensation (BEC) 
at least in dimensions greater than two. Let us consider 
this case. 



n = n o + Jdni9i/d+i(h\,z), 



(46) 



where no is the number of particles at the first quantum 
level (with k = 0) divided by V. The function gu is given 
by 



9k(\l\,z) = 



(h\z) r 



(47) 



m= 1 



BEC occurs when \j\z = 1 and then the state k = 
starts to be microscopically occupied. This gives that in 
three dimensional space we have the following condition 



.1] 



A 3 n|7| = 2.612. 



Hence, we find the critical temperature 

2nh 2 ,„UA 



k B T2 



m 



2.612 



(48) 



(49) 



This temperature is smaller for particles with fractional 
bosonic statistics than for ordinary bosons. Below T c the 
contribution to the internal energy comes from particles 
which are outside the condensate and we can easily find 
that 



U 
V 



T(d/2 + l)f(d/2+ r 



r(d/2)£(d/2) 
Hence, the specific heat is 



T 

n(k B T) ( ^ 



d/2 



(50) 



r(d/2 + l)f (d/2 + 1) fd + 2 



r(d/2)£(d/2) 
and the equation of state is 



T 

kn I j^f 



d/2 



P : 



dr(d/2 + i)e(d/2 + r 

2 r(d/2)f(d/2) 



n(k B T) 



T 

X? 



d/2 



(51) 



(52) 



The dependence on 7 is hidden in only. On eq. I|49|) 
we see that as I7I — > then — > too, and hence, the 
thermal functions l|50(l - l|52(l diverges. This expresses 
the fact that the limit 7 = is a nonanalytic boundary 
between particles with 7 > and with 7 < 0. In other 
words, this particular value of 7 establishes a true phase 
boundary between two phases of noninteracting 7-ons. 

In the next part we are going to discuss the low tem- 
perature properties of g-ons. At T = the momentum 
distribution function reduces to 



%(ei) = ~ e *)' 



(53) 



for all < g < 1. This means that for all positive g we 
have a true Fermi surface in the reciprocal space which 
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separates occupied and unoccupied states. Note also, 
that the case 3 = must be treated separately because 
then n g {ci) is infinite. Again, as in the case of 7-ons with 
7 > 0, the main contribution to the low energy properties 
is given by the particle-hole excitations around the Fermi 
surface. Performing the Sommerfeld-type expansion we 



find that (see Appendix B and also |llj) 
„ d-2 



(J- 



24 



■gCMk B TY 



(54) 



In the present case there is no linear temperature cor- 
rections to the chemical potential. We also calculate the 
internal energy 



U 



-^e g + n{e s F ) l - 2 ' d C 1 {g)(k B Tf, 
where now the Fermi energy is 

Hence, we obtain the specific heat 
C v = K g T, 



(55) 



(56) 



(57) 

where K g ~ (l/g) 2 / d ~ 1 Ci(g). The low temperature 
properties of g-ons resemble those of ideal fermions with 
modified density of states by the factor 1/g. We also see 
that the case g = can not be reached analytically from 
the case g > 0. This is due to the possibility of BEC of 
ideal bosons (g = 0). 

As an additional problem it was interesting to consider 
BEC for ideal g-ons with g < 0. As we proved rigorously 
in the Appendix C, there is no BEC in this case. So, only 
true bosons can develop a macroscopic occupation of the 
lowest (k = 0) quantum level at finite temperatures. 

Finally, we consider p-ons at low temperatures. At 
T = the distribution function (|1 T|) is 



n p (ei) = pQ(p - e;). 



(58) 



We see that, although both functions in eq. (|T7|) diverge 
at €i = the final result is well behaved; those two sin- 
gularities cancel each other. The last result means that 
p-ons possess a sharp Fermi surface with a wavevector 



(59) 




h 2 \2ripJ 



This property allows us to perform again the Sommerfeld 
(low temperature) expansion for thermodynamic quanti- 
ties. This was done in Appendix B. We note that for any 
1 < p < 00 the specific heat is linear with temperature, 



CONCLUSIONS 

In this paper we studied the thermodynamic proper- 
ties of particles obeying three different statistics. Each 
of them interpolates the distribution function between 
the BE and the FD limits. We compared both the high 
and the low temperature properties of 7-ons, g-ons and 
p-ons. We showed that 7-ons behave like fermions when 
7 > and like bosons when 7 < at all temperatures. 
However, their properties are not the same because their 
thermal functions are parameterized by 7 which changes 
the number of states. We also showed that in between 
of the boselike and the fermilike regimes, there is a true 
boundary condition (7 = 0), where the statistics is clas- 
sical. However, when we approached this limit from the 
left or the right sides the thermal functions diverge as 
I/7. In fig. 1 we draw the schematic phase diagram for 
7-ons for different values of 7. 

The thermodynamic properties of g-ons turned out to 
be quite different. At high temperatures we noted two 
classes of behavior: fermilike for g > 1/2 and boselike 
for g < 1/2. For g = 1/2 the statistics was classical. 
However, at low temperatures g-ons resemble ordinary 
fermions with modified density of states for all < g < 1 . 
This means that there must be a crossover temperature 
when the typical behavior changes its character. The 
schematic phase diagram is provided in fig. 2. 

Finally, we considered p-ons and we showed that for 
p > 1 they look like bosons at very high temperature. 
On the other hand, they resemble fermions at low tem- 
peratures if p < 00. So, again there must exist a crossover 
temperature from the bosonic to the fcrmionic behavior. 
The phase diagram for p-ons is shown in fig. 3. 

In conclusions, we showed that different routes of in- 
terpolating between FD and BE limits corresponds to 
different ways of approaching the statistical mechanics 
of particles with a nonstandard statistics. Our methods 
2 and 3 are based on the combinatorial approach, which 
is not directly related to the symmetry properties of the 
many-body wave function with respect to the particle 
transposition. In particular, the knowledge of the second 
virial coefficient is not enough to determine the full ther- 
modynamics. For example, although we know the virial 
coefficient for anyons, it is impossible distinguish which 
statistics describes them at low temperatures, 7-ons or 
g-ons. We think that none of them is the proper one. 
However, anyons in the strong magnetic field seem to be 
described by g-on statistical mechanics 



Cv = K P T, 



(60) 



where now K p ~ p-^jP 2 ' ■ So, we conclude that at 
very low temperatures p-ons resemble free fermions with 
modified density of states. 



APPENDIX A 



In this Appendix we show how to solve eq. I|10|) using 
the Lagrange inversion theorem. Suppose, that we have 



7 



a transcendental equation 

y = a + x<p{y), (61) 

and we want to find y as an unknown variable. The 
variables x and a are supposed to be given. The solution 
can be represented by infinite power series 

n JlfX — 1 

n=l ' y 

Noting, that eq. IjlOfl has the same formal structure as 
(|6*T)l with a = we use eq. (|6*2*)l to obtain eq. l(2Tjl. 

APPENDIX B 

In this Appendix we present the low temperature ex- 
pansion for any thermodynamic function calculated with 



where the constant a must be chosen properly to the 
statistics. For example, for 7-ons we should take a = 
1, g-on statistics needs a = 1/g whereas for p-ons a = 
p. At low temperatures we can approximate the second 
integral extending the lower limit of integration to —00. 
Then, changing the variable in this integral y — > — y and 
expanding both functions H([j,±kBTy) around the point 
(i we obtain 

I = a H(e)de + V C l {a)- [ '^^(k B Ty+\ (65) 
Jo *! 

where the coefficients Ci(a) are defined as follows 

Ci(a) = / dyf(n(y) + (-lY(n(-y) - a)). (66) 
Jo 

If a = 1 and n(y) is 7-ons momentum distribution 
function we obtain that the even coefficients are zero 
C2n(l) = and the odd coefficients are given by 

C 2n _i(l) = (-1)™- I (2 7 r) 2 "^(l - 2 1 - 2 "), (67) 

2n 

where Bm are the Brilouin numbers. So, 7-ons have 
almost the same expansion as ideal fermions. 

Similar results we obtain for p-ons. Using the distri- 
bution function (|17|) we find that 

Can-iGO - (-!)»- W^(l - (p+ l) 1 " 2 ") (68) 



arbitrary statistics. The only limitation is that this ther- 
mal function must be well behaved in the vicinity of 
chemical potential. Let us consider the following inte- 
gral 



1=1 H(e)n{e- /j,)de (63) 
Jo 



Introducing the new variable y = (3e — f3[i and manipulate 
with the integrals we obtain 



(64) 

I 

and Cm = 0. Note that for p — 00 we obtain the co- 
efficients for the Sommerfeld expansion of ideal bosons 
provided that they do not condense. 

For g-ons it is impossible to find the analytic formula 
for Ck(g)- We only know that C (g) = [lj and the 
rest ones are finite. 

The formalism developed above is useful in calculating 
the chemical potential and the internal energy. Namely, 
the equation for /j is 

n = — = Q / dee^-ye), (69) 
V Jo 

and for the energy 

771 /■CO 

y=^J o dee d / 2 n{t), (70) 

where 

I / jyi \ d/2 

Q = (27r) d / 2 r(d/2) \¥J ■ 

Expanding eq. (|69|l to first nontrivial order as in eq. I|65|) 

we have 

n ~ fl(a f dee 4 ' 2 - 1 + C 1 (a)(d/2 - l)(/j) d / 2 - 2 (fc s T) 2 ), 
Jo 

(72) 



J 



I = - / d V H{\ + n)n(y) + - / d V H{\ + M )(n(y) - a) + - / dyH{\ + M ) a , 

P Jo P P J-Pfi P P J-Pfi P 
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and then supposing that the correction to the chemical 
potential is small fx = e F + Sfj, we can find that 



/i = 4 1 - (d/2 - 1) 



Ci(a) /& B T 



(73) 



Similarly, we find the low temperature internal energy 
E 



V 



+ n(e a F ) d / 2 - 1 C 1 (a)(k B T) 2 



(74) 



where e a — §q^§ (e^) d ^ 2+1 . The last equation allows us to 
calculate the specific heat which is 



C v = K a T, 



(75) 



where K a = 2Q{^ l ) 1 - 2 / d C 1 {a)k 



no is the population of the lowest energy level and we 
treat it separately in case of BEC. Eq. (|78|l is absolutely 
converged for all z in the range 



(t) + l) n+1 
0<z<z R = w 



(80) 



(77 + 2) f '+ 2 ' 

On the other hand, momentum distribution function can 
be written in the form where w(£) obeys the equa- 
tion 



w v (l + w) 



l+V 



1 



,/3e 



(81) 



Condition for BEC is w = r\. Substituting this into eq. 
(|76|l and taking the maximal value of possible z = zr we 
obtain the allowed e 



APPENDIX C 



fie = 2(l + ?7)ln(l + 77) - (2 + 77) ln(2 + 77) 77 In 77. (82) 



In this Appendix we extend our discussion of <?-ons to 
the case with negative g (i.e. we suppose that g < 0). 
In principle, the value of g can be arbitrary and it seems 
very interesting to investigate how the particles look like 
with " attractive" statistical interaction. The most excit- 
ing question deals with the BEC in such system. How- 
ever, as we show below, the condensation does not occur 
in such gas. 

The equation for the momentum distribution function 
is following (cf. JTHJ) ) 

n g (e^ - (1 + (1 + 77)77 9 (e i )) 1+ '"(l + nn^))^, (76) 

where we define g — —77, with 77 > 0, and £ = e^ 6 * - 
Using the method developed in ref. 01 ( see a ls° Ap- 
pendix A) we find the solution of eq. (|76|l in the form of 
the infinite power series, i.e. 



(77177 + V + m )» 



m=0 



rn ! 



(77) 



This series is absolutely converged for £ > ^f^yjrr- 
Next, we find the equation for the total particle num- 
ber 



N 1 1 . . 



(78) 



where the function <?fc(?7, z) is given by the following series 
(77m + 77 + m) rn z m + l 



m— 



(m + l) k 



(79) 



We see that in order to have BEC e must be negative but 
we have supposed that e > 0, so we obtain contradiction 
and this proves the absence of BEC for g-ons. 
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FIGURE CAPTIONS 



Fig.l. Schematic phase diagram for 7-ons. 
Fig.2. Schematic phase diagram for p-ons. 
Fig.3. Schematic phase diagram for p-ons. 



